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Abstract. We present the first non-LTE atmosphere models for WR stars that incorporate a self-consistent solution of the 
hydrodynamic equations. The models account for iron-group line-blanketing and clumping, and compute the hydrodynamic 
structure of a radiatively driven wind consistently with the non-LTE radiation transport in the co-moving frame. We construct 
a self-consistent wind model that reproduces all observed properties of an early-type WC star (WC5). We find that the WR- 
type mass-loss is initiated at high optical depth by the so-called 'Hot Iron Bump' opacities (Feix-xvi). The acceleration of the 
outer wind regions is performed by iron-group ions of lower excitation in combination with C and O. Consequently, the wind 
structure shows two acceleration regions, one close to the hydrostatic wind base in the optically thick part of the atmosphere, 
and another farther out in the wind. In addition to the radiative acceleration, the 'Iron Bump' opacities are responsible for an 
intense heating of deep atmospheric layers. We find that the observed narrow O vi-emissions in the optical spectra of WC stars 
originate from this region. By their dependence on the clumping factor we gain important information about the location where 
the density inhomogeneities in WR-winds start to develop. 
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1. Introduction 

The modeling of the atmospheres of Wolf-Rayet (WR) stars 
with their strong stellar winds underwent a strong develop- 
ment during the last two decades. Starting with pure-he lium 
atmospheres dHilliedl 1 987albt lHamann & Schmut dfl987l) . the 
main goal was to augmen t the models by more complex atoms 
like C , N, and O (e.g. iHillierl Il989t iKoesterke & Hamannl 
1 19951) . until a major step towards a realistic determination 
of stellar parameters and emergent flux distributions was 
performed by the inclusion of iron-group lin e-blanketing 
dHilUer & Millerll998lfl999tlGrafener et alJ2002l) and clump- 
ing fsee lHam^mi&Koe^tfiAeTl998lK 

For the inclusion of line-blanketing major revisions of the 
model codes were necessary. In line-blanketed models the 
whole spectral range is covered by millions of overlapping 
lines. Due to the frequency coupling in expanding atmospheres, 
no sampling or re-ordering techniques can be appli ed (cf. th e 
technique of opacity distribution functions, e. g. ICarbonl 1 979l) . 
i.e., the complete radiation field must be calculated on a fine 
frequency grid, with each line transition at its proper place. 
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As a by-product of these numerically expensive calcula- 
tions also the radiative acceleration on the wind material is ob- 
tained. This is of special interest for WR stars, because it is 
not clear whether their strong winds can be driven by radia- 
tion pressure alone. In contrast to the thin winds of OB stars, 
WR winds show wind momenta Mvoo far above the single scat- 
tering limit L/c. This means that, on average, each photon leav- 
ing the star must be scattered several times to accelerate the 
wind to the observed terminal velocities. 

The first results from line-blanketed atmosphere calcula- 
tions have shown that the acceleration of the outer part of 
a WR wind is easily explained by radiative driving, whereas 
the models failed compl etely in the deeper l ayers where 
the mass-loss is initiated jGrafener et alJ l2000l iHerald et alJ 
200 it Idrafener et alJl2002t lHillierll2003i On t he oth er hand, 
Nu gis & Lamersl d2002l) and lLamers & Nug isl d2002h demon- 
strated by a critical-point analysis that the mass-loss rates of 
WR stars can probably be maintained, if the critical point of 
the equation of motion is located in the optically thick part of 
the atmosphere at a temperature close to 160kK. In this case 
the mass-loss is initiated by the so-called 'Hot Iron Bump' - an 
increase of the Rosseland mean opacity, which is mainly due to 
the excitation of the Fe M-shell ions (Feix-xvi). Due to their 
method of analysis, however, it remained unclear whether the 
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assumed conditions are actually met in WR atmospheres, and if 
the wind can be further accelerated beyond the optically thick 
regime. 

Due to their high wind densities and their hard radiation 
field, WR stars develop extended atmospheres in extreme 
non-LTE that require sophisticated modeling techniques. 
Their complex ionization structure is strong ly affected by 
the in terplay between different elements (see IGrafener et alJ 
120021) . It is therefore necessary to employ a full non-LTE 
treatment and a radiation transport that correctly accounts 
for line-overlaps. Moreover, WR winds show high velocities 
at large optical depths, i.e., complete atmosphere models 
including an accurate solution of the temperature structure are 
required. The numerical cost of our present models is therefore 
muc h higher than for already existent stellar wind models 



tinuum radiation' that causes problems in case of line over- 
laps. A correct ion of the continuum r adiation field is therefore 
necessary (see IPauldrach et aflll994h . In addition, the first or- 
der Sobolev approximation breaks down in the acceleration re- 
gion near the transsonic wind base - even for thin stellar winds 
JSellmaier et al Jl 99/SHOwocki & Pul Jl 999h . 

The hydrodyna mic structure of O sta r winds was obtained 
for the first time by ICastor et all Jl975l later on CAK), who 
employed a parametrization of the radiative acceleration by 
means of the so-called force multi plier function. This me thod 
has been significantly improved bv lPauldrach et alJ Jl986l) and 
describes successfully the observed properties of OB star winds 
(except for dense win ds where multiple scatter i ng effects 
become important, see lLamers & Leifhererlll993l IPuls et ail 
fl995lVinketal]l200nl) . However, due to the parametrization 
by only two parameters, the line force is subject to an averag- 
ing process, i.e., the hydrodynamic equations are not strictly 
solved. Again, for the highly structured winds of WR stars, this 
which are mostly focusing on the thin ap proach is probabl y not sufficient. The Monte-Carlo models 
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Krticka & Kubatll2C 
winds of OB stars. 

An additional important aspect for the driving of 
WR winds, which is not included in any of the previous wind 
models, is the role of density inhomogeneities. For WR stars, 
these are directly observed as time-dependent variations of 
the line profiles JLep ne & Moffat! 1 19991: Lepine et all 1 19991 
120001 iKoesterke et al.ll200*I ~ They are included in our mod- 
els in the form of a clumping factor D, which lowers the 
empirically derived mass-l oss rates by a factor Vd (see 
lHamann & Koesterkelfl998l) . i.e., it reduces the wind momen- 
tum problem for WR stars. Typically, values in the range of 
D =4-10 are obtained by a fit of the electron scattering wings 
of strong emission lines. An effect of similar importance is the 
influence of clumping on the radiative acceleration o f the wind 
materi al. T his effect was disc ussed for the first time bv lSchmutzl 
lll995|) andlSchmu tzl lll997|) and is also detected in our mod- 
els ( IGrafener & Hamannl2003l) . The clumping factor therefore 
plays a key role for the dynamics of WR winds. 

For the modeling of OB star winds, two main approaches 
are presently followed. In the first one, all atomic popu- 
lations are calculated in non-LTE (i.e., the statistical equi- 
librium equations are exactly solved) and the radiatio n field 
is treated in Sobolev approximation JPauldrach etalJl200ll 
IKrticka & Kubatll2004l) . In the second, a Monte-Carlo tech- 
nique is employed for the radiation transfer which particularly 
accounts for multiple scattering effects due to line-overlaps. 
The atomic populations of the iron-group element s, however, 
are tr eated in a modified nebular approximation JVink et al] 
Il999l) . As mentioned above, both methods are probably not 
suitable for the modeling of WR winds, because the spectrum 
formation in WR stars is dominated by non-LTE effects and 
line-overlaps. Despite of the approximations concerning the 
atomic populations, all spectral lines are treated as pure scat- 
tering opacities in the Monte-Carlo approach, i.e., they do not 
allow for photon leakage from one part of the spectrum to an- 
other (see lOwocki & Gavlevlll9^LSirr]l2004 This situatio n 
may be improved by the new technique of lLucvl (1200211200^ . 
The Sobolev approach, on the other hand, introduces an artifi- 
cial distinction between 'line radiation' and 'background con- 



approacn is probably not sufficient, lne Monte-Carlo models 
bv lvinket all i 19991) are based on an even less accurate method 
where a consistent glo bal energy budget is prov ided for a pre- 
scribed velocity field. IKrticka & Ku Mtl J2004I) have recently 
improved this situation by a full local solution of the equation 
of motion and the energy equation, i.e., their models provide 
a consistent velocity and temperature structure. Nevertheless, 
due to the neglect of multiple scatterings, these models are lim- 
ited to weaker winds. 

In the present work we overcome the problems explained 
above by combining full non-LTE atmosphere models with the 
equations of hydrodynamics. By iteration between the non- 
LTE radiation transport and the hydrodynamics we derive M, 
v(r), T(r), and the non-LTE population numbers consistently 
with the radiation field in the co-moving frame (CMF), i.e., 
we solve the radiation transport, the energy equation, and the 
equation of motion simultaneously. Moreover, clumping is ac- 
counted for. The resultant models describe the conditions in 
WR atmospheres correctly, and provide synthetic spectra, i.e., 
they allow for a direct comparison with observational material. 
Nevertheless, simplifying assumptions are necessary also in 
our models. These concern especially the assumption of a con- 
stant Doppler broadening velocity throughout the atmosphere 
and the omission of opacities, partly due to the neglect of trace 
elements like Ne, Ar, S, or P and partly due to incompleteness 
of the available data (for ions above ionization stage ix, no fine- 
structure data, and no atomic data for elements with Z > 26 are 
available). 

As already mention ed above, the critical -point analysis by 
iNugis & LamerJ J2002I) revealed that the mass-loss of early- 
type WR stars is probably driven by the 'Iron Opacity Bump' 
which is mainly due to line opacities of the Fe M-shell ions 
(Feix-xvi). By inclusion of these ions in o ur models we are 
able to check whether the assumptions of INugis & Lamersl 
(critical point at high optical depth and high temperature) are 
correct, and most importantly, if the high wind acceleration can 
be maintained in the outer part of the wind. 

In the Sect.|2]of the present work we review our standard 
atmosphere models, followed by a series of test calculations to 
clarify the role of the 'Iron Bump' opacities and the clumping 
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factor (Sect.01. Then we describe the combination of our at- 
mosphere models with the hydrodynamic equations (Sect.0}, 
and present a self-consistent hydrodynamic atmosphere model 
that reproduces all observed characteristics of an early-type 
WC star. Finally, we summarize our conclusions in Sect.[6] 

2. Model atmospheres 

In the present section we give a short overview of our 
standard atmosphere models with prescribed wind struc- 
ture, as they are routinely u s ed for the spectral analysis 
of WR stars fcoesterke et all Il992t 
Leuenha gen & H amannl Il994c 
1995 : [ Leuenhagen etaD 1 199 ' 
19981 IKoesterke et all 1200: ' 



Hamann et all 1 19921 



Koesterke & Hamann 



Hamann & Koesterke 
IGrafener et all 120021 
Hamann & Grafened l2003l) . In Sect.lTTl we describe the 
applied numerical methods, in Sect. l2.2l we review the relevant 
model parameters, and in Sect. 12.31 we summarize the atomic 
data utilized for the calculations. 

2.1. Standard models 

For our standard models we assume a spherically symmetric, 
stationary outflow with a prescribed density and velocity struc- 
ture. The model code computes the radiation field, the atomic 
populations, and the temperature structure in the expanding at- 
mosphere. The complete solution comprises three parts, which 
are iterated until consistency is obtained: The radiation trans- 
port, the statistical equilibrium equations, and the energy equa- 
tion. 

The radiation transport for the spherically expanding at- 
mosphere is formulated in the co-moving frame of refer- 
ence (CMF), neglecting aberration and advection terms (see 
iMihalas et al1ll976l) . For a fast solution and a consistent treat- 
ment of the scatter ing terms we employ th e method of variable 
Eddington factors dAuer&Mihalasll970l) . This means that the 
moment equations are solved to obtain the angle-averaged ra- 
diation field, and the numerically expe nsive ray-by-ray trans- 
fer is only calculated from time to time ( IKoesterke et all2002l 
IGrafener et aHl2002l) . The fast numerical solution allows a de- 
tailed treatment of millions of spectral lines on a fine frequency 
grid. 

The atomic populations and the electron density are deter- 
mined from the equations of statistical equilibrium. This sys- 
tem of equations is solved in line with the radiation transport 
by applicati on of the ALI formal ism (accelerated lambda it- 
eration, see lHamannlll985l Il986l) . Complex model atoms of 
He, C, O, Si, and the iron-group are accounted for. For the 
inclusion of millions of iron-group tran sitions we take advan - 
tage of the concept of super-levels (see IGrafener et aHl2002h . 
Furthermore, density inhomogeneities are accounted for by the 
assumption of small-scale clumps with a clumping factor D 
(defined as the inverse of the volume filling factor fy, see 
IHamann & Koesterkell998l) . 

The temperature structure is obtained from the assumption 
of radiative equilibrium. In the present work this constraint 
equation is decoupled from the equations of statistical equilib- 
rium. It is solved by a temperature correction procedure which 



is ba sed on the Unsold-Lucy method JUnsoldl 1195.4 iLucvl 
and has been generalized f or the application in non-LT E 
models with spherical expansion JHamann & GrafeneJ2003l) . 



2.2. Model parameters 

The model atmospheres are specified by the luminosity and ra- 
dius of the stellar core, the chemical composition of the enve- 
lope, and its density and velocity structure. The basic parame- 
ters are: The stellar core radius at Rosseland optical depth 
tr = 20, the stellar temperature (related to the luminosity 
L+ and the core radius 7?* via the Stefan Boltzmann law), the 
chemical composition, the clumping factor D = 1 //V, and, if 
the atmosphere structure is not calculated self-consistently, the 
mass-loss rate M and the terminal velocity In the latter case, 
a yS-type velocity-law of the form u(r) = Voo ( 1 - Rolrf is as- 
sumed, where Rq denotes the transition point to the hydrostatic 
domain and /? is set to one. 

At the start of the model calculation, when the density and 
velocity structure of the atmosphere are fixed, only the con- 
tinuum opacities for LTE populations with a grey tempera- 
ture structure are calculated. The value of tr = 20, which 
specifies the inner boundary of our models, therefore denotes 
a 'LTE-grey-continuum Rosseland optical depth' of 20. The 
'true' Rosseland optical depth, tr oss , of a converged non-LTE 
model with all line opacities included is usually much larger. In 
our present models (WC stars with solar Fe-abundance) values 
of Tross = 50 are typically reached at the inner boundary. 

Together with the density and velocity structure also the 
clumping factor D is prescribed at the start of the model cal- 
culation. For realistic wind models it is however necessary to 
allow for a radial variation of D. Especially for our intended 
hydrodynamic treatment a smooth wind-structure is favorable 
at large optical depth around and below the sonic point. The 
radial dependence is therefore specified in the following form. 
In the deep layers, where tr is greater than a specified value 
Ti, a smooth wind with D = 1 is assumed. Then the clumping 
factor is continuously increased until a fixed maximum value is 
reached at a specified optical depth T2. For tr < Tn, D is kept 
constant. In the present work, the given values for D refer to its 
maximum value with t\ set to 0.7 and to to 0.35. 

The parameters R*, M, Uco, and D are connected by 
the transformed radius R t oc R*(vco/ VdM) 2 ^ 3 . Models with 
the same R t show almost identical line equivalent widths 
JSchmutz et al.lll989l) . This invariance holds because the line 
emis sion in WR stars is dominate d by recombination processes 
(see IHamann & Koesterke! 19981) . 

2.3. Atomic data 

The chemical composition of the model atmospheres is given 
by mass fractions Xn e , Xq, Xq, X^ and Xp e of helium, car- 
bon, oxygen, s ilicon and iron-group elements. Compared to 
IGrafener et all J2002I) the atomic data are extended by the ions 
Fe x-xvn and by detailed photoionization cross sections for C 
and O that include dielectronic resonances. Two types of model 
atoms, with and without the extended iron model, are applied in 
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17 
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16 
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18 


30122 
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Fe v 


19 
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32 
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18 
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Cm 


40 


Fe vii 


16 
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54 


Fe vm 


10 
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Fe ix 


11 
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On 


3 


Fe x 
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761 


Om 


33 


Fe xi 
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1143 


Oiv 


25 


Fe xii 


9 


894 


Ov 


36 


Fexm 
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797 


Ovi 


15 


Fexiv 


9 


589 


Ovn 


1 


Fe xv 


7 


253 


Sim 


10 


Fe xvi 


6 


31 


Si iv 


7 


Fe xvii 


1 


1 


Siv 


1 









Table 1. Summary of the extended model atom. The iron-group 
ions (Fe) are described by a relatively small number of super- 
levels, each representing a large number of true atomic energy 
levels (sub-levels). In the present work we also apply a smaller 
model atom, where the iron-group ions are restricted to Fe n- 
Fe x with Fe x represented by one auxiliary level. 



the present work. The model atoms are summarized in Tabled 
They contain ionization stages He i-He m, C i-C v, O n-O vn, 
Si m— Si v, and Fe n-Fe x or Fe n-Fe xvn. Except for He i and 
Si in, the lowest and highest ionization stages are restricted to 
few auxiliary levels. 

The atomic level & line data are compiled from the 
following sources . Oscillator strengths are t aken from the 
Opacity Pro j ect JHummer & Mihalasl 1 19881 ISeatorJ 1 19871 
Seatonetal.l Il992t ICunto & Mendozal Il992t ISeatonetall 
19941 iThe Opacity Project Teaml Il995t). level en ergies from 
Kurucz's CD-ROM No. 23 jKurucz & Belllll995l) . iron group 
data (io nization stage s II-IX ) from Kurucz CD-ROM No. 20- 
22 (see iKuruczlf 1 99 ll \200A and the iron M-shell data (Fex- 
Fexvn) from the Opacity Project. The iron-group elements 
from Sc to Ni are comprised into one generic atom. This model 
atom is of course dominated by Fe due its high relative abun- 
dance. The extension of the generic ion from the stages Fe x to 
Fe xvii is restricted to pure Fe because the Opacity Project does 
not provide data for elements with Z > 26. 

The emission line spectra of WR stars are chiefly formed in 
line emission cascades, which follow recombination to excited 
energy levels. A realistic treatment of photoionization and re- 
combination processes is therefore mandatory. Especially for 
the ionization of carbon and oxygen ions C i-C m and O i-O v, 
the excitation of auto-ionizing, dielectronic resonances close 
to the ionization edge provides an important reaction cannel. 
In our previous work, photoionization cross sections were ob- 
tained by low-order fits to Opacity Project data, and the di- 
ele ctronic transitions w e re trea ted separately by the approach 
of iMihalas & Hummerl dl973l) . This distinction between di- 
electronic and 'pure' photoionization processes is, however, 
highly questionable. The resonances actually form complex 
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Table 2. Model parameters for our different test models 
(A,B,C, and D) and the hydrodynamic model (Hydro) com- 
pared to the values derived for WR 1 1 1 by JGrafener et alJ 
120021) . Stellar core parameters: Stellar luminosity L*, radius 
R+, temperature T*, and mass M+ (relevant only for the hy- 
drodynamic model). Wind parameters: Mass-loss rate M with 
the corresponding transformed radius R t , clumping factor D, 
and terminal wind velocity u^. Atmospheric abundances: Mass 
fractions of helium, carbon, oxygen, silicon and iron X He , Xc, 
Xo, Xsi, and Xp e . Wind efficiency: Resultant work-ratio Q = 
Wwind/^wind according to Eqs. Q and (0, and wind momen- 
tum efficiency 77 = L^/cMv^. 



features on the total cross section, which cannot be eliminated 
by a fitting process. In the present work we avoid these prob- 
lems by taking the ioniz ation cross s ection s for C and O by 
iNahar & Pradhanl dl997h and iNaharl dl999h as they are, i.e., 
they are only smoothed with a Doppler profile. Effects like 
the successive saturation of the dielectronic features when a 
continuum becomes optically thick, are therefore modeled in a 
realistic manner. The numerical cost of this treatment is, how- 
ever, rather high because the continuum opacities must be pro- 
vided on a fine frequency grid to resolve the dielectronic fea- 
tures. Further disadvantages may be due to the slightly incor- 
rect wavelengths in the data source. Moreover, the cross sec- 
tions are only available in a form where the transitions from 
one lower level to several states of the upper ion are combined. 
Photoionizations therefore always lead into the ground state of 
the upper ion. 

3. Test calculations 

As a preparation for our intended hydrodynamic modeling, we 
investigate in the present section the effect of the additional Fe 
opacities on the radiative acceleration, and determine the range 
of stellar parameters that actually allows to maintain a WR-type 
wind. For this purpose we employ models with a prescribed 
yS-type velocity law and similar parameters as the WC 5 pro- 
totype WR 1 1 1, which has a lready been analyzed in detail by 
iHillier & Millerl i 19991) and iGrafener et ai] d2002h . Although 
the modeling of WC stars is a rather complex task, these ob- 
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Fig. 1. Test models A and B (T* = 85 kK) with prescribed ve- 
locity field: Wind acceleration in units of the local gravity g{r) 
vs. atomic density as depth index. The mechanical + gravita- 
tional acceleration as implied by the adopted /Maw (black) is 
compared to the wind acceleration due to radiative + gas pres- 
sure. Whereas Model A (grey) comprises the iron group ions 
Fe n-Fe x, Model B (grey dashed) is extended by the full iron 
M-shell up to Fe xvn. 



v/km s~ 



1500 1000 300100 30 10 



1 



T" 



T" 



1 — I T 

- Model C 

- Model D 




10 12 
log 10 (n tot /crn 3 ) 



14 



16 



Fig. 2. 'Hot' test models C and D (r* = 140 kK) with pre- 
scribed velocity field: Wind acceleration as in Fig.^but for 
models with = 140 kK. For Model D (grey, dashed) the 
clumping factor has been increased from D — 10 — > 50 
and the mass-loss rate has been correspondingly reduced by 
V50/10. For this model a globally consistent energy budget 
with Q = 1 .05 is obtained whereas Model C only provides a 
fraction of Q — 0.71 of the wind energy. 



jects offer a good test case for our hydrodynamic models. Due 
to their almost homogeneous chemical composition their mass 
can be reliab ly derived from the theoretical mass-luminosity re- 
lation e.g. bv lLangerlll 19891) . i.e., the number of free parameters 
is reduced by one. 

In Sect. l3.ll we identify the line opacities of the Fe M-shell 
ions as the main contributors to the wind acceleration in deep 
atmospheric layers, and show that the WR-type mass-loss can 
only be maintained for rather high stellar core temperatures 
(around = 140 kK). Furthermore, we point out the influence 
of the clumping factor on the radiative acceleration. In Sect. 13. 21 
we demonstrate the reaction of the optical O vi-features on the 
radial dependence of the clumping factor, and show that they 
are a possible diagnostic tool for the location where the wind 
clumping sets in. 

3.1. The effect of the Fe M-shell opacities 

To investigate the effect of the additional Fe opacities we start 
with a similar model (i.e. with similar paramet ers and atomic 
data) as in our previous analysis of WR 111 jGrafener et alJ 
I2002I) . The model parameters are adapted in a way to facilitate 
the intended hydrodynamic treatment (slightly higher C and O 
abundances) and to allow for an easy comparison to the final 
model in Sect.|5] (i>«, — > 2020 kms~'). Moreover, the mass- 
loss rate is adjusted to reproduce the observed Cra/Civ line 
strengths. Note that the effective wind density is enhanced in 
our present models due to the improved treatment of dielec- 
tronic recombination (see Sect. l2.3> . In Table|2] the parameters 
of this first test model (Model A) are listed together with the 
values previously derived for WR 111 (GKH02). 

The radiative acceleration as obtained from Model A is 
shown in Fig.^ It is compared to the wind acceleration accord- 



ing to t he prescribed /?-type velocity law. In analogy to previous 
work s dG rafener et al.l200ollHerald et alJ200lllGrafener et alJ 
l2002l lHillierll2003l) . the obtained radiative force exceeds the 
force needed to accelerate the wind in the outer layers, but fails 
at large optical depth where the mass-loss is initiated. For the 
work ratio Q between the work performed by the radiation field 
and the gas expansion on the wind material 



W, 



wind 



= m f 

Jr. 



and the prescribed mechanical wind luminosity 
1 



-'wind 



M 



2 M+G 



(i) 



(2) 



we obtain a value of Q = Wwind/^wind = 0.47, i.e. the model 
provides about half of the global energy budget of the wind. 

The reason for the lack of force in the acceleration region 
becomes obvious when the atomic populations of Model A 
are inspected. The top panel of Fig.[3]shows the Fe-ionization 
structure of Model A (black). In the inner part of the atmo- 
sphere, Fe is ionized up to the highest ionization stage consid- 
ered in the model atom, i.e., in this region only opacities from 
subordinate levels contribute to the radiative force. The lack of 
force may therefore be attributed to a lack of sufficiently high 
ionization stages considered in the model atom. For C and O 
the situation is similar, however, these elements are complete 
up to their He-like ions. Fe is only considered up to the sec- 
ond M-shell ion, Fe x. Because the ionization potentials of the 
M-shell ions are increasing only in small steps from 236 eV for 
Fe ix to 489 eV for Fe xvi before a big step of 1266 eV follows 
for Fe xvii, one may expect that at least a few more M-shell 
ions could be excited if they were present in the model atom. 
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Fig. 3. Test calculations: Ground-state 
populations of the included iron-group 
ions vs. atomic density as depth index 
for models A, B, and C. Top panel: The 
'cool' models A (black) and B (grey) 
show ionization stages up to Fexiv. The 
model atom for Model A which only in- 
cludes ions up to Fex is clearly insuffi- 
cient. For both models, the M-shell ions 
(Fe rx-Fe xvn) are strongly concentrated 
towards the inner boundary. Bottom panel: 
The 'hot' Model C shows M-shell ions al- 
ready at lower densities. Ionization stages 
up to Fe xvi are excited. 



That this is indeed the case is indicated by the grey curves 
in Fig. [3] They show the populations of Model B which has 
the same parameters as Model A but an extended model atom, 
including all M-shell ions up to Fexvn (see TableQ. In this 
model, ions up to Fe xiv are excited. However, an inspection of 
Fig.[2shows that the radiative force due to the additional ions 
does not solve our problem with the wind acceleration. a ra( j is 
only increased in the hydrostatic layers, where it now consid- 
erably exceeds the gravitational attraction, i.e., no hydrostatic 
equilibrium is possible anymore. The deficit in the acceleration 
region remains. 

The apparent solution of this problem is an increase of T+, 
to push the temperature regime where the M-shell ions are ex- 
cited farther out into the acceleration region of the wind. For 
a typical WC model, the emergent spectrum is only marginally 
changed by this operation because the continuum becomes op- 
tically thick already at high wind velocities, i.e., the location of 
the photosphere is determined by the velocity s tructure of the 
wind rather than the size of the stellar core (see lHamann et al.1 
120031) . In the deeper layers however, where our problem with 
the radiative acceleration occurs, the temperature scales with 

and the ionization structure may be changed considerably. 

In the lower panel of Fig.|3]we show the Fe-populations of 
Model C where the stellar radius R+ has been correspondingly 
decreased so that = 140 kK. It is clearly visible that the bulk 
of the M-shell ions is now shifted towards much lower densi- 
ties, and ions up to Fe xvi are exited. As demonstrated in Fig.[2] 
the radiative acceleration on these ions works exactly in the 
intended way, and fills the previous gap in the acceleration re- 
gion. Moreover, the work ratio Q of this model has improved to 
a value of 0.71. Nevertheless, below and directly above the ap- 
parent opacity bump, the radiative acceleration still lies slightly 
below the gravitational attraction. Whereas the deficit in the 
deeper layers, below log(« tot /cm _1 ) = 13.5, would be bal- 
anced by the gas pressure in a hydrodynamic model, the lack of 
force around log(« tot /cm _1 ) = 12, would presumably remain. 
The latter, however, is located so far outside in the wind that 



changes in the recombination rates due to clumping can con- 
siderably affect the wind dynamics. 

This is demonstrated by Model D in Fig.|2l where the 
clumping factor has been increased from D = 10 to a value 
of 50. The mass-loss rate has been correspondingly reduced by 
a factor of V50/10 in order to maintain the transformed radius 
R t and thereby the strength of the emergent emission line ra- 
diation (see Sect. l2.2> . For this model a work ratio of 1.05 is 
obtained, i.e., it has a nearly consistent global energy budget. 
With respect to Model C, the radiative acceleration is partially 
enhanced by up to 0.3 dex. The influence of the clumping fac- 
tor and its exact dependence on the radius is therefore crucial 
for the dynamics of WR winds. 

The lack of force in Model C is presumably due to the omis- 
sion of trace ele ments like Si, P , S, Ne, and especially Ar in our 
calculations (see lHillieill2003b . i.e., the high clumping factors 
in Model D are possibly only needed as a compensation for 
this deficit. Nevertheless, we use Model D as the start model 
for our hydrodynamic modeling in Sect.|5J and admit that we 
might overestimate the clumping factor in order to achieve a 
consistent solution. 

Interestingly, with an electron temperature of 140-200 kK 
in the acceleration regi on, our Model D basical ly matches the 
conditions claimed by iNugis & Lamersl J2002I) for the initi- 
ation of an optically thick wind. However, our models indi- 
cate that the acceleration of the wind above the sonic point is 
only possible for very hi gh core temperatur es of * 140 kK. 
This is in contradiction to lNugis & Lamersl who assumed much 
larg er radii, in accordance w ith current spectral analyses (see 
also lLamers & Nugisll 2002l). These typically provide values of 
7\ < 90 kK dKoesterke & Hamannll99llGrafener et al.ll998l 
120021 ICrowfher et alJkOojT Our present value, on the other 
hand, is in accordance with stellar structure calculation s which 
also s how high core temperatures (120-130kK, see lLangeJ 
ll989llSchaerei1ll996l) . This may imply that, due to the obscu- 
ration of the stellar core by the wind opacity, spectral analyses 
tend to overestimate the stellar radius, i.e., the large derived 
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Fig. 4. Model spectra around O vi 38 1 1/34 for different clump- 
ing structures. For Model D our standard approach is used (see 
text). For Model D2 clumping is initiated farther out in the wind 
whereas for Model Dl a constant clumping factor is assumed. 
The Ovi emissions react strongly on the clumping structure. 
The farther out clumping is initiated the stronger becomes O vi. 
For Model D2 also the continuum-level is affected. 



As demonstrated in Fig.|4] Ovi3811 reacts strongly on the 
performed changes of the clumping factor. In our test, Model D 
shows slightly stronger emissions than observed. Model Dl, 
which is fully clumped, provides much weaker lines, whereas 
for Model D2, where the clumping sets in farther out, not only 
O vi is enhanced, but also the flux level of the continuum is 
affected. Although the strength of the Ovi lines will depend 
on the final model parameters (see next section), the present 
test shows that O vi is a good diagnostic tool for the clumping 
structure in deep atmospheric layers. 

4. The wind hydrodynamics 

In this section we describe the combination of our model atmo- 
spheres with the equations of hydrodynamics. The resultant at- 
mosphere models provide a self-consistent wind structure, i.e., 
for given stellar core parameters the non-LTE population num- 
bers and the atmosphere structure are derived consistently with 
the radiation field in the co-moving frame of reference (CMF). 
The solution is obtained by iteration between the non-LTE ra- 
diation transport and the hydrodynamics. Thereby, the numeri- 
cal effort is considerably enhanced with respect to the standard 
models. In Sect. l4.ll we discuss the parametrization of the ra- 
diative acceleration, which forms the link between the radiation 
transport and the hydrodynamics. The solution of the hydrody- 
namic equations is explained in Sect. 14.21 



radii do actually reflect a very extended wind-structure rather 
than a large core size. 

3.2. The radial dependence of the clumping factor 



As explained in Sect. 12.21 the radial dependence of the clump- 
ing factor is fixed at the start of our model calculations. For the 
models in Sects. |3~T1 and IBl the clumping factor is assumed to 
increase from a homogeneous wind (D = 1) for optical depths 
larger than t\ = 0.7 to a specified value (D = 10 or 50) for 
optical depths smaller than T2 = 0.35. As demonstrated in the 
previous section, clumping has a strong impact on the radia- 
tive acceleration, i.e., the results of our hydrodynamic model- 
ing will depend on the detailed choice of t\ and T2. 

For a justification of our present choice we have calculated 
test models similar to Model D but with different values of t\ 
and T2- For Model Dl the clumping factor is assumed to be 
constant throughout the atmosphere, whereas for Model D2, T] 
and T2 have been decreased (t! = 0.3, T2 = 0.1) so that the 
onset of clumping is pushed farther out in the wind. 

In the emergent spectrum these changes especially show 
up in the nar row Ovi-doublet at 3 811/34A. In previ- 
ous works (e.g iHillier & Milled 1 1999t Idrafener et aDl2002l 
ICrowther et al.l I2OO2I) these lines could not be reproduced, 
which is remarkable because of the relatively simple structure 
of the Ovi ion. That we are now able to reproduce Ovi fea- 
tures of the observed strength is due to the combination of the 
high stellar temperature and the additional blanketing by the Fe 
M-shell opacities. O vi consequently helps to trace regions rela- 
tively deep inside the envelope at temperatures around lOOkK. 



4.1. The radiative acceleration 

From the solution of the moment equations in the radiation 
transport, we obtain the Eddington flux H v on a fine frequency 
grid in the CMF. In combination with the opacity k v , which is 
given on the same grid, it is thus possible to derive the radiative 
acceleration on the wind material by direct integration 



flrad 



1 An f 00 

K V H V 

P c J 



dv. 



(3) 



Due to Doppler shifts, however, the radiative acceleration 
on spectral lines depends strongly on the velocity gradient v'. 
The radiative-acoustic wave modes arising from this depen- 
dence are chan ging the characteristic properties of the gas flow 
dAbbotdll98(j . For the application in our hydrodynamics it is 
therefore necessary to find a suitable parametrization of a ra( j 
which accounts for this effect. 

H fl975l) the acceleration 



In the CAK approach llCastor et 
on spectral lines is written in the form 



C 1 1 dv 
r 



C /ckA° 
drj ' 



(4) 



where we have used our own definition dr = pdr. In this 
formulation the parameter a reflects the distribution of line 
strengths in the relevant wavelength range, and is usually found 
to be in the range 0.6-0.7. In contrast to this, the radiative ac- 
celeration in our models, as calculated from Eq. Q, cannot be 
described by two parameters. To obtain consistent wind models 
we have to allow for depth dependent parameters C and a. By 
this, however, the solution topology due to the non-linearity in 
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v' becomes much more complicated. To avoid numerical prob- 
lems we therefore utilize a linearized form of Eq. @ in our 
calculations. 

Let us first assume that aij nes is given for a fixed velocity 
field vo(r). Then a first order expansion of Eq. (|4} with respect 
to v' gives 



n , C ( dv ° 
(1 - a) 7 d7 



+ a- 

r L 



C ( dv$ 
d7 



(5) 



Here the first term oc l/r 2 behaves as expected for optically 
thin lines (note that vq is fixed), and the second term oc j^/pr 2 
reflects the behavior of optically thick lines due to Doppler 
shifts, i.e., the parameter a indicates the fracti on of the line 
force which is d ue to optically thick lines (see Ida vlevll 199.4 
IPuls et all2OO0l) . 

For a realistic representation in the hydrodynamics, a ra d 
from Eq. (0 must be disentangled into its different contribu- 
tions. As a first step, the line acceleration aii nes must be sep- 
arated from the acceleration on continua a cont . Then the ratio 
between optically thin and optically thick lines must be deter- 
mined, so that a ra d can be written as a cont + a t hin + flthkk- Qthin 
and fltMck correspond to the first and second term in Eq. (J5J, i.e., 
they must be parametrized in a form where a t Mn « 1 /r 2 and 
tfthick 37 lr 2 - Because the acceleration on continua is also in- 
dependent from the velocity gradient, a cont can be treated anal- 
ogously to fl t hi n . Finally, a ra( j can be written in the form 



Orad - (acont + flthin) + «tl 



MGF(t) k(r)(dv\ 

-^ 2 ~ + — £• () 



where F(r) includes the contributions from continua and opti- 
cally thin lines, and Ic(t) describes the contribution from opti- 
cally thick lines. 

The first step, the separation of a cont and aii nes , is easily done 
by dividing the total opacity k v under the integral in Eq. Q 
into the parts due to continua and lines only. Then the depen- 
dence of aiines on v' is investigated to determine a tn i n and a t hick- 
For this purpose the velocity field v(r) in the radiation trans- 
port is varied, i.e., the radiation field is calculated twice for the 
same density structure and the same atomic populations but for 
different velocity distributions u{r). In practice, a variation of 
10% is applied to v(r) and v'{r). The variation of v' then causes 
changes in the derived value of aii nes and the total radiative flux. 
After normalizing the result with respect to the flux, we obtain 
a measure for the reaction of aii nes on the variation of v' alone. 
Now the fraction of the line force due to optically thick lines 
(which is equivalent to the force multiplier parameter a) can be 
calculated directly from Eq. (0 for each depth index. 

The physical meaning of a derived on this way is differ- 
ent from the standard approach, where a is interpreted as a 
characteristic material property due to the distribution of line 
strengths. The value derived here is an effective value which 
reflects the reaction of aii nes on v' including all important ra- 
diative effects like multiple scattering, line overlaps, ionization 
effects, the exact angle-dependence of the radiation field, or the 
shortening of the photon mean free path at large optical depth, 
which also ha ve partly been introduced in the standard CAK 
appro ach (see lAbbotl|[l982l IPauldrach et alJll986l IPuls et alJ 
2000). 



The parametrization of a la( j in Eq. (|SJl is fundamentally dif- 
ferent from the CAK approach (Eq.0Ji which is non-linear in 
v'. Our linearization helps to avoid numerical problems due 
to multiple solutions. However, in the CAK approach the non- 
linearity is prerequisite for the self-consistent calculation of the 
mass-loss rate. For our present calculations, we therefore need 
an additional constraint. Because F and k are continuously up- 
dated, the overall nonlinear behavior of the line force is never- 
theless retained during the iterative solution. 

4.2. Solution of the hydrodynamic equations 

The self-consistent hydrodynamic atmosphere models are ob- 
tained by iteration between the non-LTE radiation transport as 
described in Sect. 12.11 and the solution of the hydrodynamic 
equations. The radiation transport provides the radiative ac- 
celeration according to Eq. l|6} and the temperature structure, 
namely T(r), k{r), and T(t). Based on this, a hydrodynamic so- 
lution for the density and velocity structure, p(r) and v(r), is 
obtained, p and v are then fed back into the radiation transport, 
and the procedure is repeated until convergence. 

As already indicated above, T(t) as well as F(t) and £(r) 
are fixed on the T-grid, where r is given by our simple defi- 
nition dr = pdr from Eq. @. By this, the energy equation is 
implicitly included in the hydrodynamics (at large optical depth 
T is coupled to Tr oss ). This is particularly important for the op- 
tically thick winds of WR stars where a violation of radiative 
equilibrium results in extreme deviations from flux constancy, 
which would spoil the whole iteration process. 

To maintain radiative equilibrium it is further necessary that 
the total optical depth of the atmosphere remains unchanged 
when a new hydrodynamic structure is calculated. As explained 
in the following, this additional constraint is exactly what we 
need to close the system of hydrodynamic equations, i.e., by the 
inclusion of the energy equation we gain the additional con- 
straint which is required due to our linearization of a ra( j with 
respect to v' . We incorporate the conservation of optical depth 
into the hydrodynamics by choosing r as the integration vari- 
able. Due to our simple definition of r, however, we neglect the 
the actual p-dependence of flux-mean opacity. Nevertheless, 
this approach turned out to be numerically stable and efficient 
enough for our purposes. 

With a la d from Eq. l|6), the equation of motion for a station- 
ary, spherically expanding gas flow reads 



dv 
'dr 



MG(l -F) 



1 dp k 1 1 dv 



p dr 



p dr 



(7) 



where M denotes the stellar mass, G the gravitational constant, 
p the gas density, and p the gas pressure. Formulating p in 
terms of the isothermal sound speed a = ^p/p, replacing dr 
by dr/p, and applying the equation of continuity 



M = 4npvr , 

we obtain the equation of motion in the form 



(8) 



do 
d7 



_a_ ] M 



-MG(l-D + 2a 2 r-*L^, (9) 
4jtv dr 
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where T(t), 1c(t), and a(r) = tJRT(t)/h are taken from the 
preceding non-LTE model calculation. This equation must be 
solved in line with r(r) which is given by our definition of t 

dr = 1 = W 
At p M 

For given boundary values Db and r\» and a given mass-loss 
rate M, Eqs. (|9} and dlOi can be integrated directly over r. The 
outer boundary of our model atmosphere is located at a pre- 
scribed radius r\, = R max . At this point t can be arbitrarily set 
to Tb = 0. The inner boundary value t* at the stellar core radius 
R+ is taken from the previous non-LTE solution. By these two 
conditions, namely r(r*) = R+ and r(0) = R msix , we ensure that 
the optical depth of our atmosphere model is conserved. A third 
condition is obtained at the critical point of the equation of mo- 
tion, where the term in square brackets in Eq. (|9j becomes zero. 
In analogy to the CAK approach, this point marks the position 
where the flow spee d reaches the value of the fast radiative- 
acoustic wave mode jAbbotJl980l) . At this point also the right 
hand side of Eq. (|9} must be zero for dv/dr to be defined. 

After all we have three conditions to obtain two boundary 
values «b and r^, and the mass-loss rate M. For the starting 
point of the numerical integration we choose the outer bound- 
ary, which is located at radius R mdx with r = 0. By this, the 
first of our three conditions (r(0) = ^ max ) is fulfilled, and the 
values v(Q) = V&, and M remain to be determined. As start val- 
ues for and M we take the values from the previous model 
iteration. Then the numerical integration of Eqs. (|9} and dlOi 
is performed, until the critical point is reached. If the critical 
condition is not fulfilled, M is adjusted and the integration is 
repeated until this is the case. Now the numerical integration 
can be continued through the critical point to the inner bound- 
ary at t = r*, where our last condition r(r*) = must be 
fulfilled. This is achieved by an adjustment of [»«,, and repeti- 
tion of the previous procedure until r(r + ) = R+. 

Instead of using a pure numerical algorithm we employ rel- 
atively simple rules for the adjustment of u M and M. At the 
critical point, the new value for M can be read off directly from 
Eq. (|9j. «oo, on the other hand, can be adjusted by the well- 
known result from the CAK theory that the terminal velocity 
of a radiatively driven wind scales with the escape velocity at 
the stellar surface, i.e., if r{r+) + R+ we scale the terminal ve- 
locity by v esc (R*)/v esc (r). 

The numerical integration is performed by a simple Runge- 
Kutta scheme with adaptive step-size control. Here especially 
the calculation of dv/dr from Eq. is affected by numeri- 
cal difficulties. The resolution of the model grid in our non- 
LTE models (typically 70 depth points) is much lower than 
what is needed for the integration of the hydrodynamic equa- 
tions. For the extraction of T, k, and a it is thus necessary 
to employ an interpolation algorithm which provides smooth 
functions with well-de fined derivative s. We use the monotonic 
spline-interpolation by Steffen C 1990) for this purpose. Further 
difficulties occur when the critical point is approached, where 
Eq. (|9}, which has the form dv/dr = f/g, has a singularity with 
/ — » and g — > 0. Close to this point it is necessary to apply 
the rule of de l'Hopital for the calculation of dv/dr. 
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Fig. 5. Wind acceleration as in Fig.^ but for the self-consistent 
hydrodynamic model. The grey line indicates the acceleration 
flwind due to radiation + gas pressure as calculated within the 
radiation transport. a„i n d is in precise agreement with the me- 
chanical + gravitational acceleration (a mec h + «grav, black dots) 
due to the hydrodynamic velocity structure and the stellar mass. 
In addition, the radiative acceleration on free electrons (axh), on 
continua (a C ont), and the total radiative acceleration on continua 
and lines (a ra d) are indicated. 

In practice, the wind structure is iterated until changes in 
the mass-loss rate AM/M < 10 -4 are obtained, which is typi- 
cally the case after 15-30 iterations between non-LTE radiation 
transfer and hydrodynamics. The demands on the numerical 
accuracy of the non-LTE solution are extremely high. Already 
small deviations from radiative equilibrium can have fatal con- 
sequences, because the critical point reacts very sensitively on 
changes of the radiative flux. An accuracy better than 0.01% is 
therefore mandatory for the flux-level. Because of this, the non- 
LTE solution alone requires a computation time of the order of 
one day on a current personal computer, and a hydrodynamic 
model like the one presented in Sect.[5]takes a pure calculation 
time of the order of one month. 

5. Self-consistent modeling of a WCstar wind 

In this section we apply the hydrodynamic treatment described 
in Sect. 14.21 and present a self-consistent atmosphere model for 
an early-type WC star. The 'input' parameters for the hydrody- 
namic model are the same as for our test model D in Sect.|3] 
Moreover, we assume a stellar mass of 1 3.63 M R according to 
the mass-luminosity relation bv lLangell Jl989l) . The clumping 
structure is prescribed in the same way as for our previous cal- 
culations, i.e., D(tr) is related to the LTE-continuum optical 
depth tr (see Sect. l2.2t which is now updated after each itera- 
tion. 

The converged model nicely reproduces the observed pa- 
rameters of a typical WCstar (see Table|2] 'Hydro'). The ob- 
tained mass-loss rate is only slightly higher than the value 
adopted for Model D, and the resultant terminal velocity is 
very close to the canonical value of 2000 km/s. In Figs.Qand 
[8] the synthetic spectrum is compared to the observed spec- 
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Fig. 7. Hydrodynamic model: Comparison of the synthetic spectrum (thick line, grey) to the observed spectrum of WR 111 (thin 
line). Prominent spectral lines are identified. After correction for interstellar extinction, the observed flux and the model flux are 
both divided by the model continuum. Additionally, a correction for interstellar Ly a absorption is applied. The model parameters 
are compiled in Table|2](Hydro). 



trum of WR 111. Although the fit quality would not be suf- 
ficient for a detailed analysis, the model clearly displays the 
spectral characteristics of an early-type WC star. Moreover, the 
successful reproduction of Ovi 38 11/38 supports our assump- 
tions concerning the radial dependence of the clumping factor 
(see previous section), although the electron scattering wings of 
C iv 5805 indicate that the assumed maximum value of D = 50 
is probably too high. 

Further discrepancies are mainly due to the modified model 
parameters in compar ison to our previous analysis of WR 111 
jGrafener et ai1l2002l) . Due to the high assumed carbon abun- 
dance, many carbon lines appear too strong. Moreover, the 
emergent flux distribution is slightly changed, which results 
in a smaller required correction for interstellar extinction 
(Eb-v = 0.25 mag compared to the previous value of Eb-v - 
0.325 mag). In the present work, the model flux must there- 



fore be scaled down to match the observed flux distribution 
of WR 1 1 1 . A comparison of the emergent flux distributions 
is shown in Fig.|8] where the model flux has been reduced by 
0.18 dex and a distance of d — 1.58 kpc (M - m = 1 1.0 mag) is 
adopted for WR 111. 

The wind acceleration of our self-consistent model is 
shown in Fig. [5] A precise agreement between the mechanical 
+ gravitational acceleration (black dots) and the acceleration 
due to radiation and gas pressure (grey, in the background) is 
obtained. Note that the wind shows two acceleration regions, 
one in the optically thick part below a velocity of 900km/s, 
and one in the outer part of the wind. These two zones reflect 
the influence of the two well-known opacity bumps, which are 
responsible for a variety of phenomena in stell ar physics (see 
iRoeers & IriesiaJl998tlNugis & Lamersl2002l) . In the present 
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Fig. 8. Hydrodynamic model: Emergent flux distribution (grey) compared to the observed flux of WR 111 (dashed black line) 
including optical and infrared photometry (labeled blocks). A distance modulus of 11.0 mag (d = 1.58 kpc) and interstellar 
reddening (with Eb-v = 0.25 mag) are applied to the calculated spectrum (grey) and the pure continuum calculation (solid black 
line). Moreover, the model spectrum is scaled down by 0.18 dex to match the absolute value of the observed flux. 
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Fig. 6. Velocity structure of the hydrodynamic model (black) 
compared to di fferent 0-tyye velocity l aws and a double-/? law 
as proposed bv lHillier & Milled dl999h . 



case, both opacity bumps are located beyond the sonic point 
and thus contribute to the acceleration of the wind. 

In Fig. [6] the resulting velocity structure is compared to 
/?-type velocity la ws with 0-1 and = 5, and a double-/? 
law as proposed bv lHillier & Millerl dl999h (wco = 2000 kms , 
Uext = 400 kms 1 , Pi = 1, and 02 = 50). Clearly, none of the 
/Haws is able to describe the overall structure of the hydrody- 
namic model exactly. The strong acceleration in the inner part 
is best described by /3 = 1 w hereas the outer part tend s towards 
higher /?, in accordance with lLepine & MoffatN 19991) who de- 
duced high values (fi > 3) for the outer wind region from vari- 
ations in emission line-profiles. The best match, especially for 
the asymptotic b ehavior at large radii, is indeed provided by 
Iffilfier & MilierF s double-/? law. 



For the parametrization of the line force in the hydrody- 
namics, the effective force multiplier parameter a is calculated 
within our radiation transport (see Eq.|4}. In the standard pic- 
ture, a denotes the fraction of the line force which is due to 
optically thick lines, i.e., the fraction which reacts linearly on 
changes of v'. As demonstrated in Fig. [9] we obtain values for 
our WR model of -0.1 < a < 0.1 throughout the whole atmo- 
sphere. This implies that the dependence of the line force on the 
velocity gradient is very small and sometimes even shows the 
opposite behavior as expected (i.e. the line force decreases for 
increasing v'). The displayed Ostar model, in contrast, shows 
the expected behavior with a ~ 0.7 in the outer wind regions, 
and a — » at large optical depth. The abnormal values do pre- 
sumably not reflect the original distribution of line strengths. 
Instead, the effective value of a is probably influenced by com- 
plex radiative processes like line emission cascades due to re- 
combination in the outer wind regions or extreme line overlaps. 
Hence our results indicate the total breakdown of the CAK the- 
ory for WR winds, even in the outer wind regions. 

Concerning the driving mechanism of the wind, our models 
are in qualitative agreement with the ideas of lNugis & Lamersl 
d2002l) . The critical point of the equation of motion (at u c = 
32 kms -1 ) is indeed located very close to the sonic point (at 
v s = 28 kms -1 ) in the optically thick part of the envelope, be- 
low the iron opacity peak where the Rosseland mean opacity 
is increasing with radius. The radiative acceleration on the iron 
opacities overr uns the gravitational attr action close to the sonic 
radius (see also lNugis & Lamersl2002l Sect. 4.1), i.e., the star 
crosses the Eddington limit at this point and thereby initiates 
the mass loss. 

The detailed structure of our present model, however, 
shows significant differences to the results of lNugis & Lamersl 
for WR 111. T his is basically due to the high value 
iNugis & Lamersl adopted for the stellar radius (correspond- 
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ing to T+ — 8 1 kK, see lLamers Because 
iNugis & Lamersl onlv analyze the sonic point, and assume that 
the wind is radiatively driven in the outer regions, they ob- 
tain solutions for a relatively wide range of T*. Our test cal- 
culations, in contrast, indicate that for 'cool' around 85 kK 
the radiative force is not sufficient to maintain the WR-type 
mass-loss beyond the sonic point (see Model B in Sect.[3jl. For 
our 'hot' models with = 140 kK this situation improves. 
Consequently, as the temperature gradient increases propor- 
tional to r*, we obtain higher temperatures at smaller optical 
depths. The sonic point in our hydrodynamic model is therefore 
located at a relatively small flux-mean optical depth of r s = 5.4 
and a high temperature o f T s =199 kK, compared to the mod- 
els by INugis & Lamersl where r s = 15-25 and T s = 140- 
190 kK. 

In our hydrodynamic models the mass-loss rate is chiefly 
determined by the temperature-dependence of the Rosseland 
mean opacity. An increase of the wind density heats up the 
deeper layers and pushes the opacity bump, and therefore also 
the sonic point, towards lower densities. Consequently, the 
mass-loss rate is decreased, i.e. the wind is stabilized. This reg- 
ulation mechanism is most effective in parameter ranges where 
the ionization structure is sensitive to changes in the mass-loss. 
This is probably the reason why the WR models in spectral ana- 
lyses tend to lie in regions where the ionization is just changing 
(i.e. the models are most difficult to calculate). 

Another important question is how our result s fit to the pre- 
dicted propertie s of stellar structure models (e.g. iLangerl 19891 
ISchaererlll996l). A n inspection of Figs. 4 and 5 in the work 
of lSchaererl ([1996) shows that our sonic-point temperature of 
r s = 199 kK is already reached relatively deep inside the star at 
around 0.8 R*. Moreover, this point is located within the outer 
convection zone of the star, which is also caused by the iron 
opacity bump in these models. This implies that the convection 
zone in stellar structure calculations may possibly be identified 
with the acceleration zone of the WR wind, which then would 
be located at a rather small radius. However, for exact predic- 
tions stellar structure models are needed that take realistic wind 
models as an outer boundary condition, and allow for convec- 
tion under dynamical conditions. 

6. Conclusions 

In the present work we have constructed the first hydrodynam- 
ically consistent model for a WR-type (WC5) stellar wind, 
driven by radiation pressure. For this purpose we have com- 
bined our non-LTE atmosphere models with the hydrodynamic 
equations. On this way we obtain the first models for radia- 
tively driven winds that incorporate a full non-LTE radiation 
transfer in the CMF, i.e., they allow for a realistic treatment 
of all important radiative effects throughout the whole atmo- 
sphere, including large optical depths. Moreover, the models 
account for the important effects of clumping. 

By inclusion of the Fe M-shell ions (Feix-Fexvi) in our 
models, the radiative acceleration in deep atmospheric layers is 
considerably enhanced. We show that, for suitable stellar pa- 
rameters, the radiative force on these opacities is able to drive 
an optically thick wind. In particular, very high stellar tempera- 
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Fig. 9. The effective force multiplier parameter a as calcu- 
lated within the radiation transport vs. radial depth index (1: 
outer boundary, 70: inner boundary). For the hydrodynamic 
WR model (black) values close to zero are obtained, indicat- 
ing a weak depende nce of the radiative force o n v'. A model 
for the O star £ Pup jGrafener & HamannlEool shows values 
around 0.7 in the outer part of the wind, as expected from the 
CAK theory. 



tures (T* « 140 kK) are needed to maintain t he wind accelera- 
tion ab ove the sonic point. In accordance with lNugis & Lamersl 
the sonic point in our hydrodynamic model is located 
at large optical depth below the 'Hot Iron Bump' at a very 
high temperature of T s = 199 kK. The resultant wind struc- 
ture shows two acceleration regions. A steep velocity increase 
in deep atmospheric layers due to the Fe M-shell ions, and a 
shallow increase in the outer part of the wind due to iron-group 
ions of lower excitation in combination with C and O. The lat- 
ter is in accordance with the shallow velocity laws deduced by 
iLepine & Moffatl i 19991) from line-profile variations. 

Despite of the fact that the wind is optically thick at the 
sonic point, we find strong deviations from the standard CAK 
theory also in the outer wind regions. Here especially the re- 
action of the radiative force on changes of u' shows an abnor- 
mal behavior, which is presumably due to complex radiative 
processes like line emission cascades after recombination, or 
extreme line overlaps. Moreover, we find a strong dependence 
of the radiative force on wind clumping, an effect which is ne- 
glected in previous models for radiatively driven winds. 

The emergent flux, as calculated from our hydrodynamic 
model, clearly shows the spectral characteristics of an early- 
type WC star, i.e., the wind density and the terminal wind ve- 
locity, which are both obtained self-consistently from the hy- 
drodynamics, are of the correct order of magnitude. Moreover, 
we are able to reproduce the observed narrow O vi emissions 
in the optical spectrum of WC stars. These features are excited 
by a combination of heating due to the Fe M-shell ions and the 
high stellar temperature, which is necessary for our wind mod- 
eling. The O vi lines also react very sensitively on the detailed 
clumping structure, i.e., they help to trace deep atmospheric 
layers where the clumping of WC star winds sets in. 
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When our results are compared to stellar structure calcu- 
lations, the question about the location of the 'surface' of a 
WR star gains more and more in importance. Stellar structure 
models for WR stars show an outer convection zone, which 
largely determines the extension of the stellar envelope. The 
convection however, is caused by the same increase of the mean 
opacity that causes the driving of the WR wind in our models. A 
temperature around 200 kK, corresponding to our sonic -point 
temperature, is only reached deep inside the stellar envelope 
at the bottom of the outer convection zone. For a full under- 
standing of the structure and the mass-loss of WR stars models 
are therefore needed, that combine stellar structure and atmo- 
sphere calculations and allow for convection under dynamic 
conditions. 
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